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Introduction 

Exclusive semileptonic decays of heavy mesons play an important role in the determina- 
tion and over-constraining of the Cabibbo-Kobayashi-Maskawa mixing matrix. The CKM 
element Vet has been extracted [|I| from B -^ D*lV and B — > DW using heavy quark sym- 
metry p|, while the element Vub has been estimated from B -^ ttW and B —* plv rates 
using various models. In both cases, the normalization and shape of the relevant hadronic 
form factors influence the extracted value of the CKM angle. For Vcb-, the normalization of 
the B — i> D'^*'> matrix element at zero recoil is provided by heavy quark symmetry. How- 
ever, typical extrapolations to this point use ad-hoc parameterizations of form factors that 
introduce theoretical uncertainties comparable to the statistical uncertainties @,^. This is 
especially unfortunate since the uncertainty in Vd, feeds into unitarity-triangle constraints 
from CP violation observed in the kaon system as the fourth power 0. For Vub neither 
the normalization nor the shape is well known. The normalization near zero recoil may be 
obtained from lattice simulations or by combining heavy quark and chiral symmetries with 
measurements of related semileptonic decays in the charmed and bottom sector 0], but a 
parameterization away from zero recoil is necessary to compare to experimental data. 

Some progress in describing the shape of such form factors has recently been made in 
the form of model-independent parameterizations 0J^ based on QCD dispersion relations 
and analyticity P,p!0|. These dispersion relations lead to an infinite tower of upper and 
lower bounds that can be derived by using the normalizations of the form factor F{ti) at 



a fixed number of kinematic points ti as input ||10|- [1^] . When the normalization is known 
at several points (say, five or more for B -^ nlv), the upper and lower bounds are typically 
so tight they look like a single line. A natural question then arises: What is the most 
general form consistent with the constraints from QCD? The answer to this question is the 
parameterization of reference Q. For a generic form factor F{t) describing the exclusive 
semileptonic decay of a i? meson to a final state meson H as a function of momentum- 
transfer squared t = [ps — PhY-, the parameterization takes the form 

1 oo 

n«) = p(j)^ !:«.*«' . (1) 

where 0(t) is a computable function arising from perturbative QCD. The function P{t) de- 
pends only on the masses of mesons below the BH pair-production threshold that contribute 
to BH pair-production as virtual intermediate states. The variable z{t\ to) is a kinematic 
function of t defined by 



l-z{t;to) \U-to ' ^^ 

where t+ = {Mb + Mh)"^ is the pair-production threshold and to is a free parameter that 
is often P,P|-P!^ taken to be t_ = {Mb — M^Y, the maximum momentum-transfer squared 
allowed in the semileptonic decay B -^ HlV. The coefficients a^ are unknown constants 
constrained to obey 

oo 

EK)'<1 ■ (3) 

fc=0 



The kinematic function z(t;to) takes its minimal physical value Zmin at t = t_, vanishes 
at t = to, and reaches its maximum z^ax at t = 0. Thus the sum J^O'k z^ is a series 
expansion about the kinematic point t = to- For B -^ D*lV with to = t_, the maximum 
value of z is z^ax = 0.06, and the series in Eq. (|l]) can be truncated while introducing only 
a small error Q]. The value z^ax can be made even smaller by choosing an optimized value 
< to < t_ p|. In that case, most form factors describing B -^ DlV and B -^ 0*117 can be 
parameterized with only one unknown constant to an accuracy of a few percent (assuming 
the normalization at zero recoil given by heavy quark symmetry). Thus the continuous 
function F{t) has been reduced to a single constant, for example the value of the form factor 
F{t = 0) at maximum recoil. For B -^ ttIv, the maximum value of z is z^ax = 0.52, but 
even in this case Eqs. (|l|) and (^ severely constrain the relevant form factor [irii[T2| . 

This remarkable constraining power can be traced to the existence of a naturally small 
parameter Zmax that arises algebraically from a conformal map. In this paper we attempt to 
trace the physical origin of Zmax in the hope of developing some intuition about the physics 
underlying the analyticity constraints of Eqs. (|I|) and @. Further, we will incorporate two 
generalizations that lead to a significantly stronger constraint on the observable B -^ -kIv 
form factor. 



Physical Basis for a Small Parameter 

To understand heuristically why there is a small parameter associated with semileptonic 
heavy meson decays, consider for the moment a form factor F{t) in the decay B —>■ DW, 
and take to = t_. Crossing symmetry tells us the analytic continuation of the form factor 
F{t) that describes semileptonic decay for < t < t_ also describes BD pair production 
for t > t^. Fig. 1 shows the general features one expects for F{t) in the region < t < oo. 
The form factor has a cut due to pair-production beginning at t = t+, as well as a series 
of poles from bound Sc-type states in the vicinity of t+. It varies rapidly near these poles, 
then falls smoothly from its peak values near t ~ t+ to its minimum values near t ~ 0. It is 
not essential to our argument that the form factor decreases monotonically as t approaches 
zero, only that the variation in F{t) over the semileptonic region < t < t_ is determined 
by the distance to the branch cut t = t^ and the magnitude of the form factor near the 
branch point, -F(t+). For fixed F(t+), F{t) varies more slowly over the semileptonic region 
as t-/t+ decreases, while for fixed t-/t-^., F{t) also varies more slowly in the semileptonic 
region as -F(t+) decreases. Both observation and QCD perturbation theory imply that the 
rate of BD pair-production cannot be arbitrarily large for t > t_^., and combined with the 
fact that t_|_ ^ t_ for B —>■ DW, we expect the variation of F{t) over the semileptonic region 
to be small. We wish to associate the small parameter Zmax with this variation. 

Suppose the form factor can be roughly described in the physical semileptonic region by 

F(t) ~ , ^° , , (4) 

^ ^ (t+ - t)P ^ ^ 

where Fq is a constant. A reasonable measure of the variation of F[t) over the physical 
region for semileptonic decay is 




FIG. 1. 
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The magnitude of a generic form factor F{t) as a function of t. The pair-production 
and the semileptonic endpoint t_ are shown schematically. 
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(5) 



This measure depends only on the kinematic thresholds t+, t_ and the power p. For p = 1 



the form factor is pole dominated and 6f is similar to the Shifman-Voloshin parameter [|13 
(Mb — MdY/{Mb + Md)'^ ~ 1/4. However, comparison with Eq. (|^) reveals that 6f can be 
identified with Zmax only if p = 1/2 giving 



5f 






(6) 



This value of p leads to the small value of Zmax for B —>■ D. Other decays such as -B — ;► plV, 
D — »• K*lT7, etc., have larger values of Zmax, with the largest occurring for B —>■ ttITj. Even 
for this extreme case, z^ax ~ 1/2 is small enough to provide a useful expansion parameter. 
On the face of it, this value of p seems rather surprising. After all, we know bound 
states exist and will contribute to form factors like poles. On the other hand, the dispersion 
relation relies on quark-hadron duality and perturbative QCD. In perturbative QCD the 
fundamental degrees of freedom are quarks and gluons, there are no bound states at any 
finite order in perturbation theory to couple to the pair-produced fermions, and the form 
factor has no poles. Indeed, at leading order in the Parton Model, the B -^ 0*117 form 



factors have the form |T^ of Eq, 



F{t) 



(7) 



with p = 1/2. Given that there are bound states in nature, how can the perturbative 
QCD results be trustworthy? Certainly perturbative QCD cannot be used directly in the 



semileptonic region. However, the perturbative calculation of pair-production should be 



reliable as long as a large region of momentum transfer is smeared over ||T5[, or integrated 
over with smooth weighting functions. By constraining the magnitude of the form factor in 
the pair-production region, the perturbative analysis indirectly constrains the shape of the 
form factor in the semileptonic region. 

For B -^ 7iW, the kinematically allowed region t_ is much larger, and the heuristic 
discussion above applies less clearly. An explicit derivation is required to see that, even in 
this case, pair-production constraints allow F{t) in the semileptonic region to be expanded 
in powers oi z < z^ax ■ 

Moments of the Dispersion Relation 

For a general semileptonic decay B — > HlV, the heuristic discussion of the previous 
section can be made concrete by considering the two-point function for the vector or axial 
vector currents J = q'j'^b, qj^j^b that arise in the charged current decay of b-hadrons, 

n7(g) = z / Ae'''"(0|TJ^(x)J^'^(0)|0) 

= (g''g^-gV'^)nJ(g^)+^^^n^(g2) . (8) 

The polarization functions Uf (g^) do not fall off fast enough at large g^ for an unsubtracted 
dispersion relation to be finite. However, derivatives of the polarization functions do fall 
fast enough at high g^ for finite dispersion relations to exist. As we wish to constrain 
hadronic form factors by a perturbative calculation, it is useful to define the derivatives of 
the polarization functions at g^ = where the partonic amplitude is well behaved, far from 
the physical region for BH pair production. At g^ = the n^'^-derivative of the polarization 
tensor for the current J'^ is 

(n) _ 1 d^^'WJjO) 

^^ ~ 3r(n + 3) 9(g2)"+2 

1 a"+iHj(0) 1 9'^+2H^(0) 

~ 3r(n + 2) 9(g2)"+i ~ 3r(n + 3) 9(g2)„+2 



IT Jo t"+3 

where J = V,A for vector and axial vector currents respectively. This dispersion relation 
relates the computation of Xj at the unphysical value g^ = to the weighted integral 
over the pair- production region of the imaginary part of Hj(g^). The higher the moment 
n the more the integral is weighted near the pair production threshold, so we expect the 
calculation to be most reliable for low moments where the smearing is largest. 

It is straightforward to determine the Xj i^i perturbative QCD. For a ratio of quark 
masses u = mq/mb, a one- loop (leading order) calculation of the vector current correlator 
gives 



Xv W 



3(r(n + 3))' 



+ 



27r2m2"+2r(2n + 6) 
1-M 



-F(n + l,n + 3;2n + Q:l - u^] 



A{n + 2) 



n + V 
uF{n + 2, n + 4; 2n + 7; 1 - M^) - F{n + 2, n + 3; 2n + 7; 1 - m^ 



(10) 



where -F(a, b; c; ^) is a hypergeoinetric function. The same expression results for the axial 
current after the substitutio 
the expressions simplify to 



current after the substitution u —> —u, i.e. Xa ("") = Xv i~u)- For a massless quark nig = 



^^'^^^^ 47r2m^+M^ + l)(^ + 2)(n + 4) " ^^^^ 



The correlator Eq. (||) has also been computed at two loops |jT6|], i.e. 0{as). The higher 
order corrections result in a 25% increase [|l^] in x^^K^)- 



Since production of BH hadrons is a subset of total hadronic production, the perturbative 
calculation x*^"^ serves to constrain the analytically continued form factors for B ^ H 
decay. More precisely, the partonic computation provides an upper bound to the smeared 
contributions of poles and cuts above the pair production threshold. The contribution of 
poles below threshold will also influence the variation of a given form factor F{t) in the 
semileptonic region, and must be considered separately. While sub-threshold contributions 
are not a fundamental aspect of the dispersion relation approach (for example, form factors 
in D — > ttIv are analytic below the D — ir threshold), they are not accounted for by the 
perturbative calculation and must be properly handled when present [|17|. 

We now turn to relating Xj to F{t). This is accomplished by inserting a sum over 
intermediate states into Imnj(g^), 

Imn-j(g^) = I J ^^^^;fj'^'^^ ^^^^ + ■ ■ ■ , (12) 

where the sum is over polarizations of H and the ellipsis denotes strictly positive contri- 
butions from the B*, higher resonances and multi-particle states. In terms of a calculable 
kinematic function k{t) arising from two-body phase space and the Lorentz structure asso- 
ciated with the form factor F{t), we may substitute the inequality 

hmIV}{t)>k{t)\F{t)\' , (13) 

into Eq. (|^) to get the contribution to the hadronic moment Xj (hadronic) from the form 
factor F(t) of interest, 

(n)., , , rii f^ kit) |F(t)P , ,, 

xTihadromc) > -^ / ^ \IJ ^' , (14) 

where nj is the isospin degeneracy of the BH pair. We rely on perturbative QCD at the 
unphysical point g^ = (or equivalently, on global duality for suitably smeared production 
rates) to assert that hadronic and partonic expressions for x are equal. Then for each n, 

x':j'\hadronic) = X?{u) , (15) 

where Xj i"^) is the n*^ moment as computed in perturbative QCD. Therefore we have that 



(n) / \ /, i-n+3 



and hence 



-/ rft|/i(")(t)F(t)|2<l , (17) 



vr Jt+ 



where h^'^\t) = J-^^jn) t (3+")/2^ The argument of the square root is positive since the 

integrand came from a production rate. 

The inequahty of Eq. (|17D makes clear how the perturbative calculation constrains the 
magnitude of the form factor in the pair-production region. To constrain the form factor in 
the semileptonic region < t < t_, we would like to find functions '^k{t) that are orthonormal 
with respect to the integral Eq. (p!7|), 



1 r°° 

- / dtUe 



ifkit) ip*{t) = 6kj , (18) 



and that vanish somewhere in the semileptonic region, say at < to ^ ^-- We could then 
expand h'^'^\t)F{t) in terms of these basis functions and use Eq. (p!7|) to bound the expan- 
sion coefficients []. If h^"\t)F(t) turned out to be analytic in t outside the pair production 
region, its expansion would be equally valid in the semileptonic region, and we would have a 
parameterization of F{t) in terms of unknown, but bounded, expansion coefficients. Unfor- 
tunately, neither h^'^\t) nor F{t) are in general analytic away from the pair production cut. 
The kinematic factor h^^\t) has explicit poles at t = and the form factor F{t) may also 
have poles arising from the contribution of bound states that can interpolate between the 
current J and the BH pair. For example, the experimentally accessible form factor /+(t) in 
B -^ nlV has a pole at t = Mb* coming from the contribution of the B* resonance. 

Fortunately, a simple pole at t = tp can be eliminated by multiplying by z{t; tp). Rewrit- 
ing z (as defined in Eq. (0) ) as 

4t;Q ^ ^^,1^,^^ , (19) 

makes it clear that z(t; tp) vanishes at t = tp and has magnitude one in the pair-production 
region, \z(t; tp)\ = 1 for t > t+. We can therefore construct a quantity with no poles outside 
the pair-production region by multiplying h^'^\t) and F{t) by factors of z{t;tp) for each 
pole at tp. To make an analytic function outside the pair-production region, we generally 
also need to eliminate square-root branch cuts in h^^\t) that arise from factors of the H- 



meson three-momentum by dividing by Jz{t,tJ). The elimination of poles and cuts from 



/i(")(t) by /i(")(t) -^ P{t)h^''\t), where P{t) is a product of ;z(t;0)'s and y'z(t,t_)'s, can be 
automatically accomplished by replacing 



^ We can choose the expansion coefficients ak to be real so that only TZe[ipk{t) ip*{t)] need vanish 
for k ^ j since it is this expression that arises in \J2 '^k^k\'^ • In a more general case where the 
Ofc are complex, our results go through unchanged if the inner product in Eq. ( [l^ ) is redefined as 
<f,g>= lim,^o ^ /t°? dt [fit + ie)g*{t + ie) + /(t - ie)g*{t - ie)]. 



1 -z{t; 0) _ 

l~^ t ~ (7^ 



t_ -t 



The elimination of poles from F{t) by F(t) -^ P{t)F{t) is accomplished by multiplying 
by a product P{t) = Iljz{t]tj) for each contributing sub-threshold resonance of invariant 
mass-squared tj. Since P{t) and P{t) have unit modulus along the pair production cut 
(the integration region in Eq. (p!7|)), the well-behaved quantity P(t)h'^'^\t)P{t)F{t) obeys 
the same relation, 

1 r°° ~ , -, 

- dt\P{t)h'^^\t)P{t)F{t)\^ <1 . (21) 

71 Jt+ 

Whereas P{t) may be viewed as a technical device to smooth out the kinematic function 
h^'^\t), P(t) contains essential information about the resonance structure of F{t) in the 
unphysical region t_ < t < t+. In Eq. (0), the poles in F{t) above threshold are constrained 
by the perturbative calculation; in Eq. (|2lD , the poles below threshold are accommodated 
by P{t). Both sets of poles influence the shape of F{t) in the semileptonic region. Since 
P(t) depends only on the position of the poles below threshold and not on the residues, it 
applies for arbitrarily strong or weak residues. Therefore, we should not be surprised that 
the eventual effect of a non-trivial function P{t) is to weaken the constraint on F{t). 
In terms of orthonormal functions (fk{t) satisfying Eq. (|l^), the expansion 



P(t)/.W(t)P(t)F(t) = $:arV.W (22) 

fe=0 

combines with Eq. (pT]) to yield 

oo „ 

E(4"0 <1 ' (23) 

k=0 

valid for moderate values of n > 0. As the expression given in Eq. (^2]) is valid everywhere 
outside the cut in the complex t plane the form factor F{t) in the region of semileptonic 
decay < t < t_ is 

^W = ^7:T77;;T7:T7;7:TE4"V.(t) • (24) 



p{t)h(-){t)p{t) t 



:0 



The sum is over A; > because by construction P{t)h'^'^\t)P{t)F{t) has no poles for 
t <t+ . All that remains is to find the orthogonal polynomials (pk{t)- This a math problem 
that can be accomplished by a change of variables. In the complex t plane, the integration 
contour may be viewed as a segment from +00 to t = t+ just below the cut and a segment 
from t = t+ to +CXD just above the cut. Defining y = yjt — t+ maps the line segments just 
above and below this cut onto the real y axis. The y axis in turn can be mapped onto the 
unit circle by the bilinear transformation ^(t; to) = {y — V^o — t+)/{y + a/^o ~ t+)- This is 



precisely the change of variables in Eq. (^. Since z^ = e'"^ are orthonormal functions on 
the unit circle, we can work backwards to find 

(t) = 1 ( t+-to \ '/' ( VU^-y/U^To ^ ' 

The expansion in orthonormal basis functions is simply a Taylor series in z{t] Iq)^. However, 
the variable z does little to aid the development of physical intuition, so we continue to 
work with the momentum transfer t. Contact with previous literature [^|[ can be made by 
identifying 



</)(-)(z(t; to)) = {^U^t + ^U - to) hf^j mh^'^Kt) , (26) 

choosing n = 0, setting tg = t- (or, in the case and language of reference p, tg = (1 — 
N)t+ + Nt_ ), and expressing the "Blaschke factors" pr8|-pU[ z{t;tp) composing P(t) in 
terms of z{t;t_) and z{tp]tJ). With these identifications Eq. ( p4D becomes Eq. (|l]) with 
4){z) = (j)'^^\z) and ai. = a)^' . 

Parameterizations for Semileptonic Form Factors 

We are primarily interested in constraining form factors that describe the decay of B 
mesons, although the formalism applies equally well to A^ baryons, or even D and K mesons 
if n(g^) is evaluated at an appropriate spacelike q^. Specialization to a particular decay and 
form factor requires an explicit computation of the functions. For a pseudo-scalar final 
meson H ot a. vector meson H* with polarization e, the various form factors in semileptonic 
B decay may be defined by 

{H*ip',e)\V'^\B{p)) = ige^'^^^elp'pp, 

{H*{p\t)\A>^\B{p)) = fe*'^ + {e*-p)[a^{p + pr + a.{p-pr] 

{H{p')\V^\B{p))= Uip + p'Y + f^P-p'r , (27) 



where it is useful to also define 



F, ' 



M 



H 



^iU - t){t- - t)a+ - ^(t - M| + Ml)f 



(28) 



with t = {p — p'Y- It is straightforward to determine the k{t) function associated with each 
of the form factors, 

fc.(^) = ^Q'[(^-M(i-^-)r^' • (29) 

For the form factors whose contribution to the rate is unsuppressed by the lepton mass, the 
indices s,p and w are given by 



kg{t) 


s = 5 


p = l 


w = 3 


kpAt) 


s = 4 


p = 2 


w = 1 


kf{t) 


s = 3 


p = l 


w = 1 


kuit) 


s = 4 


p = 2 


w = 3 



(30) 



The (j)^'^\t) functions defined in Eq. ( p6D for each form factor are 

1/4 



<Pt\t-M 



ni 



t. 



2^37rx!r^ V^H 



t 



tn 



t^ 



t 



(3+n+p) 



(Vi+-i+ 7t+-to)(^t+-i+ •^t+ - t_)"'/2(t+ - t) 



w/4 



(31) 



Our previous discussions allow us to see that each particular form factor Fi has the functional 
form 



Ht) 



1 



Ear^(t;to) 



for each moment n and expansion point to? where z may be expressed as 



(32) 



z{t]tQ) = 



\/t+ — t — Y^t+ — to 



(33) 



and X]('2 



{")n2 



< 1. 



The Blaschke factors Pi{t) depend on the masses of sub-threshold resonances. For B — >■ 
DW and B — * D*lV from factors, the masses of the relevant i?c-type resonances can be 
rather accurately estimated from potential models |2l],^. Using the results of reference 
2T| in Eq. (|l^), the Blaschke factors for the form factors / and Fi are 



P/(t)=P^,(t) 

= z{t; {6.730 GeVf)z{t; {6.73QGeVy)z{t; {7.135 GeVy)z{t; {7.U2GeVf) 

while for the form factors g and /_,_ they are 

Pg{t) = z{t; {Q.337GeVy)z{t; {Q.899GeVy)z{t; {7.012 GeV)'^)z{t; {7.280 GeVf) 

and 

Pf^ (t) = z{t; (6.337 GeVy)z{t; (6.899 GeVy)z{t; (7.012 GeVf) , 



(34) 



(35) 



(36) 



respectively. For B — > ttIu, there is only one resonance and the form factor /+(t) has the 
Blaschke factor P{t) = z{t; {5.325 GeVy). For D — ^ iilV, there are no resonances, and 
Pit) = 1. 

For the lowest moment {n = 0) and to = t-, Eqs. (|3TD - (|33D reproduce the results of 
references |^, while for arbitrary to they reproduce the results for mesons given in p|. For 
higher moments (n > 0) and a given form factor Fi, they imply 



10 



k=0 



ri [t) fc=o 



(0) 

Xj 



\ Xj 



yt7^+v^ 



E^r^t-M 



k=0 






(37) 



It is a simple matter to compute the coefficients c) and match powers of z. The condition 



fc=0 



(n) 



X) [o-k ) < 1 then implies 



(38) 



.(") 



for each n. For a given set of Cj one may find that the higher moments provide tighter 

constraints on the ttj than the lowest moment alone. This is the case for B ^ n, for the 
lowest few moments. However, for B ^ D the higher moments do not improve the bounds 
imposed by the lowest moment. 



Constraints for B -^ nlu 



To make use of Eq. ( ^2]) in extracting Vub from experimental measurements of i? ^ tiIu 
we should keep only a finite number of parameters a^ and compute the maximal truncation 
error from the omission of higher order terms in the series. This truncation error can be 
minimized by optimizing to and thereby decreasing Zmax ID- Applying Eq. (|38D further 
decreases the truncation error by restricting the range of higher order parameters (of course, 
it also restricts the range of the parameters we keep). Both effects decrease the number of 
parameters required to determine /+ (t) at a given level of accuracy. 

Experimental data on 5 — > ttW is not yet available to precisely describe the shape of 
f+{t) and so for now we will simply illustrate the utility of the higher moments. Imagine 
f+{t) is known at some fixed number of kinematic points. For concreteness we choose the 
lattice-inspired values BHH /+(21GeV2) = 1.7 and /+(OGeV^) = 0.5 to fix af,°^ and 



,(0) 



in terms of a. 



(0) 



Varying Og 



subject to the zeroth moment constraint J2{0'k j'^ < 1 



0.5 to fix Oq 

(0)x2 



maps out the envelope of parameterizations consistent with the n = dispersion relation. 
Varying Og subject to the constraint of the n*'^ moment results in a smaller allowed range 
for 02 , for the first few n. Since the allowed range of f+{t) is proportional to the allowed 
range of 02 ? a relative reduction in the range of 02 leads to the same relative reduction in 
the width of the envelope. 



^ This procedure is equivalent to the method of computing upper and lower bounds by forming 
determinants of inner products often used in the literature |9-12,26|. However, Eq. (R^) also 
dictates the shape of curves allowed inside the envelope. 



11 



The one-loop results for the allowed range of as are shown in the second column of 
Table 1. We have used a pole quark mass mi, = Mb — A corresponding to A = 0.4 GeV 
p5| (varying A by ±0.1 GeV results in no more than a 6% change in the bounds), and fixed 
to = t_. Note that the n = 3 bounds are tighter than the n = result by a factor of three. 
In fact, the n = 6 bound (not listed) is better by a factor of four. However, x^") receives 



known |16| corrections from two-loop perturbative graphs {O (as)) and nonperturbative 
matrix elements. We use values of the condensates < uu > |iGeV = (—0.24 GeV)^ and 
^ <G ' ^ 



' ixv 



Gf"' >= 0.02 ± 0.02 GeV* from reference 27 



IGeV 

The third and fourth columns of 



Table 1 show the perturbative and condensate correction factors PpJ^t and PcoLd defined by 



X^-\twoloop) = (1 + pgi, + p^::l,)x^-\oneloop) 



(39) 



Beyond n = 3 the sum of these corrections approaches 100% of the leading result and the 
reliability of the calculation becomes questionable. Even at n = 3, the corrections are large 
enough that one might worry about the convergence of the operator product expansion. 
We will use n = 2 in our examples. The allowed range of as using the two-loop result, 
including condensates, is shown in the fifth column of Table 1. Although weaker than the 
corresponding one-loop result, the two-loop bounds for n > remain significantly more 
constraining than the n = bounds. 



n 


a2{oneloop) 


(n) 
Ppert 


(n) 
rcond 


al{twoloop) 


af\HQS) 



1 
2 
3 


-0.77 < as < 0.79 
-0.45 < as < 0.47 
-0.33 < as < 0.34 
-0.26 < as < 0.28 


0.23 
0.32 
0.38 
0.42 


0.01 
0.06 
0.14 
0.28 


-0.78 < as < 0.79 
-0.52 < as < 0.53 
-0.44 < as < 0.45 
-0.40 < as < 0.41 


-0.54 < as < 0.58 
-0.40 < as < 0.41 
-0.36 < as < 0.36 
-0.35 < as < 0.35 



Ah 



Table 1. Values of the zeroth moment parameter as consistent with the n moment 
constraint. From left to right is the constraint on as at one loop, the relative two-loop 
perturbative and condensate contributions Pp2.f and p^J^d' the constraint on as at two 
loops, and the two-loop constraint on as including the B*ti spin symmetry contribution as 
described in the text. 

The ellipsis in Eq. (^) includes the contribution of intermediate B*7t states to the 
hadronic moments. In the heavy b-quark limit, the form factor g^, defined by 



ig^e'^'^'^^eo^PpP^ 



{n{p')\V^\B*ip,e)) 
is related to /+ by spin symmetry when the n meson is soft {i.e. near t = t. 

9*{v-p') 



-^f+(vp')\l + 0(—)] 
Mb ^ m'^ 



(40) 



(41) 



where v = p/Mb. Including the contribution of the B*tt state (similar contributions were 
used in [^] for the B ^ B elastic form factor) modifies the constraint of Eq. (^Tj) to 



1 

TC Jt 



dt 



P{t)X{t)<P^-} Uit)\' + |P(t)X(t)0W gM' 



< 1 



(42) 
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where X{t) = [(t+ - to)/(i+ - i)]^/^/(V^+ - i + Vi+-to)- One finds that 0^, = 0^, so 
expanding 

1 oo 

9*it) = ——M7TT.t''^itM' , (43) 



P(t)0^"Ht) ,t^o 



leads to the constraint 



E(4"^r + (4"^)'<i • (44) 

fc=0 

We cannot relate all of the b]^ to the a}^ because the b-quark spin symmetry is only valid 
near t = t^. One expects that the normalization and first derivative of g^ and /+ at t_ obey 
the heavy quark relation to ~ 10% for the physical B mass, so 

,(n) 2 0("Ht,) („) 
On = -T 7-^ On 



^0 



MBd^MtS" 



V2 1+ l"'\ ^1") 



M 



B , 



2 



.;-' ^ 472/^4"' + V2(l + /g) 4"' . (45) 



to ±10%. Taking 90% of the absolute values of the right-hand sides of Eq. ^ gives lower 
bounds on |6o I ^'^'^ l^i 1; that can be inserted into Eq. ( PI ) (this constraint applies as well to 
heavy-heavy systems, where it could be further improved by including the B*D* intermediate 
state). The improvement on the range of 0/2 from the inclusion of this additional hadronic 
final state is shown in the last column of Table 1 (using two- loop amplitudes). The effect is 
~ 30% for n = 0, decreasing to ~ 15% for n = 3. 

The ra = bounds on /+(t) at two loops, without the application of spin symmetry 
(column five of Table 1), are shown as the outer, dashed, pair of curves in Fig. 2, while the 
bounds from the n = 2 improvement, including the contribution from i?*7r, are given by the 
inner, solid, pair of curves. The constraints arising from the n > moments are a dramatic 
improvement over the lowest order n = constraint alone. Note that even the n = 1 bounds 
in the last column of Table 1 represent an improvement over that of n = alone by nearly 
a factor of two. 

The quickly weakening constraint near t = t_ would be absent had we included a nor- 
malization at zero recoil. Fig. 3 shows the spin-improved, n = 2 constraint assuming an 
additional normalization, f+{t_) = 6. It gives a feel for how tightly f+{t) can be constrained 
given measurements of /+(t) at only three points. In this case, /+(t) is determined to better 
than ±15% over the entire kinematic region. These plots are presented only to provide a 
feel for what is needed to describe f+(t) to a given accuracy. In practice, one should fit 
the parameterization Eq. (^) to experimental data and theoretical (lattice or heavy quark 
symmetry) predictions to extract Vub- The work presented here should improve the precision 
of such extractions significantly. 
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FIG. 2. Upper and lower bounds on the B -^ nlv form factor /+(t). The Dashed hnes 
are the bounds arising from the n = moment while the solid lines are the bounds arising from 
the n = 2 moment. The form factor has been fixed at two points by the lattice inspired values 
/-(-(21GeV ) = 1.7 and /+(0) = 0.5. The allowed regions are the interiors of the dashed or solid 
pairs of curves. 



14 




t (GeV^2) 



t_ 



FIG. 3. Bounds on the B -^ vr/z/ form factor /+(t) arising from three normahzation points, 
f+{t-) = 6.0, /+(21 GeV^) = 1.7 and /+(0) = 0.5. The constraints are from moments of n < 2 
as discussed in the text, including the B*tt hadronic intermediate state using bottom-quark spin 
symmetry at t = t-. 



15 



Conclusions 

We have discussed constraints on semileptonic form factors arising from analyticity and 
dispersion relations, and identified the physical origin of the naturally small parameter Zmax- 
The shape and magnitude of an analytically continued form factor in the pair-production 
region has significant impact upon the size and shape of the form factor in the semileptonic 
region. Perturbative constraints on pair-production therefore constrain semileptonic decay 
amplitudes. Such constraints take the form of parameterizations of form factors in terms 
of a small number of unknown, but bounded, constants a^. Only a small number of a^ are 
needed to describe a given form factor because the parameterizations arise from truncated 
expansions in a kinematic variable z{t; to) that is surprisingly small, \z{t; to)| < Zmax = 0.065 
for B -^ D*lV and \z{t\ to)| < Zmax = 0.52 for B -^ nlV. 

We have traced the smallness of Zmax to the dual nature of the hadronic and partonic 
descriptions of QCD. Whereas in the hadronic description the form factor is a sum of poles 
and less singular terms, in the partonic description the form factor is always less singular 
than a pole. Duality then implies that the variation of the form factor F{t) over the physical 
region < t < t_ for the semileptonic decay B -^ Hlv is characterized by 

F(t_) + F(0) ^ y^/M^+^/M^j ^"^'^^ • ^^^^ 

The small values of z^ax are due to the square-root dependence on meson masses, which is 
in turn a consequence of the absence of poles in the partonic description of QCD. 

We applied these ideas in a discussion of analyticity constraints in terms of the 
momentum-transfer variable t. We improved the constraints on the /+(t) form factor by 
~ 15 — 30% by including the contribution of the B* -^ ttIu form factor using b-quark spin 
symmetry to relate it to the B -^ ttW form factor at zero recoil. We also improved the 
analyticity constraints by considering higher moments of the dispersion relation. As an 
illustrative example, we fixed the normalization of the B — > irW form factor /+(t) at two 
kinematic points to lattice-inspired values, and plotted the envelope of allowed parameteri- 
zations (Fig.2). The higher moment constraints result in upper and lower bounds that are 
roughly twice as tight as the traditional, zeroth- moment result. Given a third normalization 
at zero recoil, the form factor is determined over the entire kinematic range to ±15%, as 
shown in Fig. 3. This suggests that a parameterization using one normalization (possibly 
given by lattice or heavy quark symmetry predictions) and two free parameters could even- 
tually be used for a model-independent extraction of Vub with roughly 30% theoretical errors 
arising from the t dependence. 
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